3.72 \(\int \frac{d+e x^2+f x^4}{x^2 (a+b x^2+c x^4)^2} \, dx\)

Optimal. Leaf size=399 \[ -\frac{x \left (c x^2 \left (-a b e-2 a (c d-a f)+b^2 d\right )+a \left (\frac{b^3 d}{a}+a (b f+2 c e)-b (b e+3 c d)\right )\right )}{2 a^2 \left (b^2-4 a c\right ) \left (a+b x^2+c x^4\right )}-\frac{\sqrt{c} \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x}{\sqrt{b-\sqrt{b^2-4 a c}}}\right ) \left (\frac{12 a^2 c e-a b^2 e-4 a b (a f+4 c d)+3 b^3 d}{\sqrt{b^2-4 a c}}-a b e-2 a (5 c d-a f)+3 b^2 d\right )}{2 \sqrt{2} a^2 \left (b^2-4 a c\right ) \sqrt{b-\sqrt{b^2-4 a c}}}-\frac{\sqrt{c} \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x}{\sqrt{\sqrt{b^2-4 a c}+b}}\right ) \left (-\frac{12 a^2 c e-a b^2 e-4 a b (a f+4 c d)+3 b^3 d}{\sqrt{b^2-4 a c}}-a b e-2 a (5 c d-a f)+3 b^2 d\right )}{2 \sqrt{2} a^2 \left (b^2-4 a c\right ) \sqrt{\sqrt{b^2-4 a c}+b}}-\frac{d}{a^2 x} \]

[Out]

-(d/(a^2*x)) - (x*(a*((b^3*d)/a - b*(3*c*d + b*e) + a*(2*c*e + b*f)) + c*(b^2*d - a*b*e - 2*a*(c*d - a*f))*x^2
))/(2*a^2*(b^2 - 4*a*c)*(a + b*x^2 + c*x^4)) - (Sqrt[c]*(3*b^2*d - a*b*e - 2*a*(5*c*d - a*f) + (3*b^3*d - a*b^
2*e + 12*a^2*c*e - 4*a*b*(4*c*d + a*f))/Sqrt[b^2 - 4*a*c])*ArcTan[(Sqrt[2]*Sqrt[c]*x)/Sqrt[b - Sqrt[b^2 - 4*a*
c]]])/(2*Sqrt[2]*a^2*(b^2 - 4*a*c)*Sqrt[b - Sqrt[b^2 - 4*a*c]]) - (Sqrt[c]*(3*b^2*d - a*b*e - 2*a*(5*c*d - a*f
) - (3*b^3*d - a*b^2*e + 12*a^2*c*e - 4*a*b*(4*c*d + a*f))/Sqrt[b^2 - 4*a*c])*ArcTan[(Sqrt[2]*Sqrt[c]*x)/Sqrt[
b + Sqrt[b^2 - 4*a*c]]])/(2*Sqrt[2]*a^2*(b^2 - 4*a*c)*Sqrt[b + Sqrt[b^2 - 4*a*c]])

________________________________________________________________________________________

Rubi [A]  time = 2.20293, antiderivative size = 399, normalized size of antiderivative = 1., number of steps used = 6, number of rules used = 4, integrand size = 30, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.133, Rules used = {1669, 1664, 1166, 205} \[ -\frac{x \left (c x^2 \left (-a b e-2 a (c d-a f)+b^2 d\right )+a \left (\frac{b^3 d}{a}+a (b f+2 c e)-b (b e+3 c d)\right )\right )}{2 a^2 \left (b^2-4 a c\right ) \left (a+b x^2+c x^4\right )}-\frac{\sqrt{c} \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x}{\sqrt{b-\sqrt{b^2-4 a c}}}\right ) \left (\frac{12 a^2 c e-a b^2 e-4 a b (a f+4 c d)+3 b^3 d}{\sqrt{b^2-4 a c}}-a b e-2 a (5 c d-a f)+3 b^2 d\right )}{2 \sqrt{2} a^2 \left (b^2-4 a c\right ) \sqrt{b-\sqrt{b^2-4 a c}}}-\frac{\sqrt{c} \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x}{\sqrt{\sqrt{b^2-4 a c}+b}}\right ) \left (-\frac{12 a^2 c e-a b^2 e-4 a b (a f+4 c d)+3 b^3 d}{\sqrt{b^2-4 a c}}-a b e-2 a (5 c d-a f)+3 b^2 d\right )}{2 \sqrt{2} a^2 \left (b^2-4 a c\right ) \sqrt{\sqrt{b^2-4 a c}+b}}-\frac{d}{a^2 x} \]

Antiderivative was successfully verified.

[In]

Int[(d + e*x^2 + f*x^4)/(x^2*(a + b*x^2 + c*x^4)^2),x]

[Out]

-(d/(a^2*x)) - (x*(a*((b^3*d)/a - b*(3*c*d + b*e) + a*(2*c*e + b*f)) + c*(b^2*d - a*b*e - 2*a*(c*d - a*f))*x^2
))/(2*a^2*(b^2 - 4*a*c)*(a + b*x^2 + c*x^4)) - (Sqrt[c]*(3*b^2*d - a*b*e - 2*a*(5*c*d - a*f) + (3*b^3*d - a*b^
2*e + 12*a^2*c*e - 4*a*b*(4*c*d + a*f))/Sqrt[b^2 - 4*a*c])*ArcTan[(Sqrt[2]*Sqrt[c]*x)/Sqrt[b - Sqrt[b^2 - 4*a*
c]]])/(2*Sqrt[2]*a^2*(b^2 - 4*a*c)*Sqrt[b - Sqrt[b^2 - 4*a*c]]) - (Sqrt[c]*(3*b^2*d - a*b*e - 2*a*(5*c*d - a*f
) - (3*b^3*d - a*b^2*e + 12*a^2*c*e - 4*a*b*(4*c*d + a*f))/Sqrt[b^2 - 4*a*c])*ArcTan[(Sqrt[2]*Sqrt[c]*x)/Sqrt[
b + Sqrt[b^2 - 4*a*c]]])/(2*Sqrt[2]*a^2*(b^2 - 4*a*c)*Sqrt[b + Sqrt[b^2 - 4*a*c]])

Rule 1669

Int[(Pq_)*(x_)^(m_)*((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4)^(p_), x_Symbol] :> With[{d = Coeff[PolynomialRemainde
r[x^m*Pq, a + b*x^2 + c*x^4, x], x, 0], e = Coeff[PolynomialRemainder[x^m*Pq, a + b*x^2 + c*x^4, x], x, 2]}, S
imp[(x*(a + b*x^2 + c*x^4)^(p + 1)*(a*b*e - d*(b^2 - 2*a*c) - c*(b*d - 2*a*e)*x^2))/(2*a*(p + 1)*(b^2 - 4*a*c)
), x] + Dist[1/(2*a*(p + 1)*(b^2 - 4*a*c)), Int[x^m*(a + b*x^2 + c*x^4)^(p + 1)*ExpandToSum[(2*a*(p + 1)*(b^2
- 4*a*c)*PolynomialQuotient[x^m*Pq, a + b*x^2 + c*x^4, x])/x^m + (b^2*d*(2*p + 3) - 2*a*c*d*(4*p + 5) - a*b*e)
/x^m + c*(4*p + 7)*(b*d - 2*a*e)*x^(2 - m), x], x], x]] /; FreeQ[{a, b, c}, x] && PolyQ[Pq, x^2] && GtQ[Expon[
Pq, x^2], 1] && NeQ[b^2 - 4*a*c, 0] && LtQ[p, -1] && ILtQ[m/2, 0]

Rule 1664

Int[(Pq_)*((d_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4)^(p_.), x_Symbol] :> Int[ExpandIntegrand[(d*x
)^m*Pq*(a + b*x^2 + c*x^4)^p, x], x] /; FreeQ[{a, b, c, d, m}, x] && PolyQ[Pq, x^2] && IGtQ[p, -2]

Rule 1166

Int[((d_) + (e_.)*(x_)^2)/((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[b^2 - 4*a*c, 2]}, Di
st[e/2 + (2*c*d - b*e)/(2*q), Int[1/(b/2 - q/2 + c*x^2), x], x] + Dist[e/2 - (2*c*d - b*e)/(2*q), Int[1/(b/2 +
q/2 + c*x^2), x], x]] /; FreeQ[{a, b, c, d, e}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - a*e^2, 0] && PosQ[b^
2 - 4*a*c]

Rule 205

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[a/b, 2]*ArcTan[x/Rt[a/b, 2]])/a, x] /; FreeQ[{a, b}, x]
&& PosQ[a/b]

Rubi steps

\begin{align*} \int \frac{d+e x^2+f x^4}{x^2 \left (a+b x^2+c x^4\right )^2} \, dx &=-\frac{x \left (a \left (\frac{b^3 d}{a}-b (3 c d+b e)+a (2 c e+b f)\right )+c \left (b^2 d-a b e-2 a (c d-a f)\right ) x^2\right )}{2 a^2 \left (b^2-4 a c\right ) \left (a+b x^2+c x^4\right )}-\frac{\int \frac{-2 \left (b^2-4 a c\right ) d+\frac{\left (b^3 d-a b^2 e+6 a^2 c e-a b (5 c d+a f)\right ) x^2}{a}+\frac{c \left (b^2 d-a b e-2 a (c d-a f)\right ) x^4}{a}}{x^2 \left (a+b x^2+c x^4\right )} \, dx}{2 a \left (b^2-4 a c\right )}\\ &=-\frac{x \left (a \left (\frac{b^3 d}{a}-b (3 c d+b e)+a (2 c e+b f)\right )+c \left (b^2 d-a b e-2 a (c d-a f)\right ) x^2\right )}{2 a^2 \left (b^2-4 a c\right ) \left (a+b x^2+c x^4\right )}-\frac{\int \left (\frac{2 \left (-b^2+4 a c\right ) d}{a x^2}+\frac{3 b^3 d-a b^2 e+6 a^2 c e-a b (13 c d+a f)+c \left (3 b^2 d-a b e-2 a (5 c d-a f)\right ) x^2}{a \left (a+b x^2+c x^4\right )}\right ) \, dx}{2 a \left (b^2-4 a c\right )}\\ &=-\frac{d}{a^2 x}-\frac{x \left (a \left (\frac{b^3 d}{a}-b (3 c d+b e)+a (2 c e+b f)\right )+c \left (b^2 d-a b e-2 a (c d-a f)\right ) x^2\right )}{2 a^2 \left (b^2-4 a c\right ) \left (a+b x^2+c x^4\right )}-\frac{\int \frac{3 b^3 d-a b^2 e+6 a^2 c e-a b (13 c d+a f)+c \left (3 b^2 d-a b e-2 a (5 c d-a f)\right ) x^2}{a+b x^2+c x^4} \, dx}{2 a^2 \left (b^2-4 a c\right )}\\ &=-\frac{d}{a^2 x}-\frac{x \left (a \left (\frac{b^3 d}{a}-b (3 c d+b e)+a (2 c e+b f)\right )+c \left (b^2 d-a b e-2 a (c d-a f)\right ) x^2\right )}{2 a^2 \left (b^2-4 a c\right ) \left (a+b x^2+c x^4\right )}-\frac{\left (c \left (3 b^2 d-a b e-2 a (5 c d-a f)-\frac{3 b^3 d-a b^2 e+12 a^2 c e-4 a b (4 c d+a f)}{\sqrt{b^2-4 a c}}\right )\right ) \int \frac{1}{\frac{b}{2}+\frac{1}{2} \sqrt{b^2-4 a c}+c x^2} \, dx}{4 a^2 \left (b^2-4 a c\right )}-\frac{\left (c \left (3 b^2 d-a b e-2 a (5 c d-a f)+\frac{3 b^3 d-a b^2 e+12 a^2 c e-4 a b (4 c d+a f)}{\sqrt{b^2-4 a c}}\right )\right ) \int \frac{1}{\frac{b}{2}-\frac{1}{2} \sqrt{b^2-4 a c}+c x^2} \, dx}{4 a^2 \left (b^2-4 a c\right )}\\ &=-\frac{d}{a^2 x}-\frac{x \left (a \left (\frac{b^3 d}{a}-b (3 c d+b e)+a (2 c e+b f)\right )+c \left (b^2 d-a b e-2 a (c d-a f)\right ) x^2\right )}{2 a^2 \left (b^2-4 a c\right ) \left (a+b x^2+c x^4\right )}-\frac{\sqrt{c} \left (3 b^2 d-a b e-2 a (5 c d-a f)+\frac{3 b^3 d-a b^2 e+12 a^2 c e-4 a b (4 c d+a f)}{\sqrt{b^2-4 a c}}\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x}{\sqrt{b-\sqrt{b^2-4 a c}}}\right )}{2 \sqrt{2} a^2 \left (b^2-4 a c\right ) \sqrt{b-\sqrt{b^2-4 a c}}}-\frac{\sqrt{c} \left (3 b^2 d-a b e-2 a (5 c d-a f)-\frac{3 b^3 d-a b^2 e+12 a^2 c e-4 a b (4 c d+a f)}{\sqrt{b^2-4 a c}}\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x}{\sqrt{b+\sqrt{b^2-4 a c}}}\right )}{2 \sqrt{2} a^2 \left (b^2-4 a c\right ) \sqrt{b+\sqrt{b^2-4 a c}}}\\ \end{align*}

Mathematica [A]  time = 1.47287, size = 444, normalized size = 1.11 \[ \frac{-\frac{2 x \left (b^2 \left (c d x^2-a e\right )+a b \left (a f-c \left (3 d+e x^2\right )\right )+2 a c \left (a \left (e+f x^2\right )-c d x^2\right )+b^3 d\right )}{\left (b^2-4 a c\right ) \left (a+b x^2+c x^4\right )}+\frac{\sqrt{2} \sqrt{c} \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x}{\sqrt{b-\sqrt{b^2-4 a c}}}\right ) \left (a b \left (e \sqrt{b^2-4 a c}+4 a f+16 c d\right )-2 a \left (-5 c d \sqrt{b^2-4 a c}+a f \sqrt{b^2-4 a c}+6 a c e\right )+b^2 \left (a e-3 d \sqrt{b^2-4 a c}\right )-3 b^3 d\right )}{\left (b^2-4 a c\right )^{3/2} \sqrt{b-\sqrt{b^2-4 a c}}}+\frac{\sqrt{2} \sqrt{c} \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x}{\sqrt{\sqrt{b^2-4 a c}+b}}\right ) \left (a b \left (e \sqrt{b^2-4 a c}-4 a f-16 c d\right )+2 a \left (5 c d \sqrt{b^2-4 a c}-a f \sqrt{b^2-4 a c}+6 a c e\right )-b^2 \left (3 d \sqrt{b^2-4 a c}+a e\right )+3 b^3 d\right )}{\left (b^2-4 a c\right )^{3/2} \sqrt{\sqrt{b^2-4 a c}+b}}-\frac{4 d}{x}}{4 a^2} \]

Antiderivative was successfully verified.

[In]

Integrate[(d + e*x^2 + f*x^4)/(x^2*(a + b*x^2 + c*x^4)^2),x]

[Out]

((-4*d)/x - (2*x*(b^3*d + b^2*(-(a*e) + c*d*x^2) + a*b*(a*f - c*(3*d + e*x^2)) + 2*a*c*(-(c*d*x^2) + a*(e + f*
x^2))))/((b^2 - 4*a*c)*(a + b*x^2 + c*x^4)) + (Sqrt[2]*Sqrt[c]*(-3*b^3*d + b^2*(-3*Sqrt[b^2 - 4*a*c]*d + a*e)
+ a*b*(16*c*d + Sqrt[b^2 - 4*a*c]*e + 4*a*f) - 2*a*(-5*c*Sqrt[b^2 - 4*a*c]*d + 6*a*c*e + a*Sqrt[b^2 - 4*a*c]*f
))*ArcTan[(Sqrt[2]*Sqrt[c]*x)/Sqrt[b - Sqrt[b^2 - 4*a*c]]])/((b^2 - 4*a*c)^(3/2)*Sqrt[b - Sqrt[b^2 - 4*a*c]])
+ (Sqrt[2]*Sqrt[c]*(3*b^3*d - b^2*(3*Sqrt[b^2 - 4*a*c]*d + a*e) + a*b*(-16*c*d + Sqrt[b^2 - 4*a*c]*e - 4*a*f)
+ 2*a*(5*c*Sqrt[b^2 - 4*a*c]*d + 6*a*c*e - a*Sqrt[b^2 - 4*a*c]*f))*ArcTan[(Sqrt[2]*Sqrt[c]*x)/Sqrt[b + Sqrt[b^
2 - 4*a*c]]])/((b^2 - 4*a*c)^(3/2)*Sqrt[b + Sqrt[b^2 - 4*a*c]]))/(4*a^2)

________________________________________________________________________________________

Maple [B]  time = 0.045, size = 1575, normalized size = 4. \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((f*x^4+e*x^2+d)/x^2/(c*x^4+b*x^2+a)^2,x)

[Out]

1/4/a*c/(4*a*c-b^2)/(-4*a*c+b^2)^(1/2)*2^(1/2)/(((-4*a*c+b^2)^(1/2)-b)*c)^(1/2)*arctanh(c*x*2^(1/2)/(((-4*a*c+
b^2)^(1/2)-b)*c)^(1/2))*b^2*e+4/a*c^2/(4*a*c-b^2)/(-4*a*c+b^2)^(1/2)*2^(1/2)/(((-4*a*c+b^2)^(1/2)-b)*c)^(1/2)*
arctanh(c*x*2^(1/2)/(((-4*a*c+b^2)^(1/2)-b)*c)^(1/2))*b*d+1/4/a*c/(4*a*c-b^2)/(-4*a*c+b^2)^(1/2)*2^(1/2)/((b+(
-4*a*c+b^2)^(1/2))*c)^(1/2)*arctan(c*x*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2))*b^2*e+4/a*c^2/(4*a*c-b^2)/(-4
*a*c+b^2)^(1/2)*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctan(c*x*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2))*
b*d-3/4/a^2*c/(4*a*c-b^2)/(-4*a*c+b^2)^(1/2)*2^(1/2)/(((-4*a*c+b^2)^(1/2)-b)*c)^(1/2)*arctanh(c*x*2^(1/2)/(((-
4*a*c+b^2)^(1/2)-b)*c)^(1/2))*b^3*d-3/4/a^2*c/(4*a*c-b^2)/(-4*a*c+b^2)^(1/2)*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c
)^(1/2)*arctan(c*x*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2))*b^3*d+5/2/a*c^2/(4*a*c-b^2)*2^(1/2)/(((-4*a*c+b^2
)^(1/2)-b)*c)^(1/2)*arctanh(c*x*2^(1/2)/(((-4*a*c+b^2)^(1/2)-b)*c)^(1/2))*d-5/2/a*c^2/(4*a*c-b^2)*2^(1/2)/((b+
(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctan(c*x*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2))*d-3*c^2/(4*a*c-b^2)/(-4*a*c+
b^2)^(1/2)*2^(1/2)/(((-4*a*c+b^2)^(1/2)-b)*c)^(1/2)*arctanh(c*x*2^(1/2)/(((-4*a*c+b^2)^(1/2)-b)*c)^(1/2))*e-3*
c^2/(4*a*c-b^2)/(-4*a*c+b^2)^(1/2)*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctan(c*x*2^(1/2)/((b+(-4*a*c+b^2
)^(1/2))*c)^(1/2))*e-1/2/a/(c*x^4+b*x^2+a)*c/(4*a*c-b^2)*x^3*b*e-3/2/a/(c*x^4+b*x^2+a)/(4*a*c-b^2)*x*b*c*d+1/2
/a^2/(c*x^4+b*x^2+a)*c/(4*a*c-b^2)*x^3*b^2*d+c/(c*x^4+b*x^2+a)/(4*a*c-b^2)*x^3*f+1/2/(c*x^4+b*x^2+a)/(4*a*c-b^
2)*x*b*f-1/a/(c*x^4+b*x^2+a)*c^2/(4*a*c-b^2)*x^3*d-1/2/a/(c*x^4+b*x^2+a)/(4*a*c-b^2)*x*b^2*e+1/2/a^2/(c*x^4+b*
x^2+a)/(4*a*c-b^2)*x*b^3*d-1/2*c/(4*a*c-b^2)*2^(1/2)/(((-4*a*c+b^2)^(1/2)-b)*c)^(1/2)*arctanh(c*x*2^(1/2)/(((-
4*a*c+b^2)^(1/2)-b)*c)^(1/2))*f+1/2*c/(4*a*c-b^2)*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctan(c*x*2^(1/2)/
((b+(-4*a*c+b^2)^(1/2))*c)^(1/2))*f+1/(c*x^4+b*x^2+a)/(4*a*c-b^2)*x*c*e-d/a^2/x+c/(4*a*c-b^2)/(-4*a*c+b^2)^(1/
2)*2^(1/2)/(((-4*a*c+b^2)^(1/2)-b)*c)^(1/2)*arctanh(c*x*2^(1/2)/(((-4*a*c+b^2)^(1/2)-b)*c)^(1/2))*b*f+c/(4*a*c
-b^2)/(-4*a*c+b^2)^(1/2)*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctan(c*x*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c
)^(1/2))*b*f-1/4/a*c/(4*a*c-b^2)*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctan(c*x*2^(1/2)/((b+(-4*a*c+b^2)^
(1/2))*c)^(1/2))*b*e-3/4/a^2*c/(4*a*c-b^2)*2^(1/2)/(((-4*a*c+b^2)^(1/2)-b)*c)^(1/2)*arctanh(c*x*2^(1/2)/(((-4*
a*c+b^2)^(1/2)-b)*c)^(1/2))*b^2*d+3/4/a^2*c/(4*a*c-b^2)*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctan(c*x*2^
(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2))*b^2*d+1/4/a*c/(4*a*c-b^2)*2^(1/2)/(((-4*a*c+b^2)^(1/2)-b)*c)^(1/2)*arc
tanh(c*x*2^(1/2)/(((-4*a*c+b^2)^(1/2)-b)*c)^(1/2))*b*e

________________________________________________________________________________________

Maxima [F]  time = 0., size = 0, normalized size = 0. \begin{align*} \frac{{\left (a b c e - 2 \, a^{2} c f -{\left (3 \, b^{2} c - 10 \, a c^{2}\right )} d\right )} x^{4} -{\left (a^{2} b f +{\left (3 \, b^{3} - 11 \, a b c\right )} d -{\left (a b^{2} - 2 \, a^{2} c\right )} e\right )} x^{2} - 2 \,{\left (a b^{2} - 4 \, a^{2} c\right )} d}{2 \,{\left ({\left (a^{2} b^{2} c - 4 \, a^{3} c^{2}\right )} x^{5} +{\left (a^{2} b^{3} - 4 \, a^{3} b c\right )} x^{3} +{\left (a^{3} b^{2} - 4 \, a^{4} c\right )} x\right )}} - \frac{-\int \frac{a^{2} b f +{\left (a b c e - 2 \, a^{2} c f -{\left (3 \, b^{2} c - 10 \, a c^{2}\right )} d\right )} x^{2} -{\left (3 \, b^{3} - 13 \, a b c\right )} d +{\left (a b^{2} - 6 \, a^{2} c\right )} e}{c x^{4} + b x^{2} + a}\,{d x}}{2 \,{\left (a^{2} b^{2} - 4 \, a^{3} c\right )}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((f*x^4+e*x^2+d)/x^2/(c*x^4+b*x^2+a)^2,x, algorithm="maxima")

[Out]

1/2*((a*b*c*e - 2*a^2*c*f - (3*b^2*c - 10*a*c^2)*d)*x^4 - (a^2*b*f + (3*b^3 - 11*a*b*c)*d - (a*b^2 - 2*a^2*c)*
e)*x^2 - 2*(a*b^2 - 4*a^2*c)*d)/((a^2*b^2*c - 4*a^3*c^2)*x^5 + (a^2*b^3 - 4*a^3*b*c)*x^3 + (a^3*b^2 - 4*a^4*c)
*x) - 1/2*integrate(-(a^2*b*f + (a*b*c*e - 2*a^2*c*f - (3*b^2*c - 10*a*c^2)*d)*x^2 - (3*b^3 - 13*a*b*c)*d + (a
*b^2 - 6*a^2*c)*e)/(c*x^4 + b*x^2 + a), x)/(a^2*b^2 - 4*a^3*c)

________________________________________________________________________________________

Fricas [B]  time = 79.0598, size = 28044, normalized size = 70.29 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((f*x^4+e*x^2+d)/x^2/(c*x^4+b*x^2+a)^2,x, algorithm="fricas")

[Out]

1/4*(2*(a*b*c*e - 2*a^2*c*f - (3*b^2*c - 10*a*c^2)*d)*x^4 - 2*(a^2*b*f + (3*b^3 - 11*a*b*c)*d - (a*b^2 - 2*a^2
*c)*e)*x^2 + sqrt(1/2)*((a^2*b^2*c - 4*a^3*c^2)*x^5 + (a^2*b^3 - 4*a^3*b*c)*x^3 + (a^3*b^2 - 4*a^4*c)*x)*sqrt(
-((9*b^7 - 105*a*b^5*c + 385*a^2*b^3*c^2 - 420*a^3*b*c^3)*d^2 - 2*(3*a*b^6 - 40*a^2*b^4*c + 150*a^3*b^2*c^2 -
120*a^4*c^3)*d*e + (a^2*b^5 - 15*a^3*b^3*c + 60*a^4*b*c^2)*e^2 + (a^4*b^3 + 12*a^5*b*c)*f^2 - 2*((3*a^2*b^5 -
13*a^3*b^3*c - 12*a^4*b*c^2)*d - (a^3*b^4 - 6*a^4*b^2*c - 24*a^5*c^2)*e)*f + (a^5*b^6 - 12*a^6*b^4*c + 48*a^7*
b^2*c^2 - 64*a^8*c^3)*sqrt((a^8*f^4 + (81*b^8 - 918*a*b^6*c + 3051*a^2*b^4*c^2 - 2550*a^3*b^2*c^3 + 625*a^4*c^
4)*d^4 - 4*(27*a*b^7 - 351*a^2*b^5*c + 1197*a^3*b^3*c^2 - 550*a^4*b*c^3)*d^3*e + 6*(9*a^2*b^6 - 132*a^3*b^4*c
+ 484*a^4*b^2*c^2 - 75*a^5*c^3)*d^2*e^2 - 4*(3*a^3*b^5 - 49*a^4*b^3*c + 198*a^5*b*c^2)*d*e^3 + (a^4*b^4 - 18*a
^5*b^2*c + 81*a^6*c^2)*e^4 + 4*(a^7*b*e - (3*a^6*b^2 + 5*a^7*c)*d)*f^3 + 6*((9*a^4*b^4 + 3*a^5*b^2*c + 25*a^6*
c^2)*d^2 - 2*(3*a^5*b^3 - 4*a^6*b*c)*d*e + (a^6*b^2 - 3*a^7*c)*e^2)*f^2 - 4*((27*a^2*b^6 - 108*a^3*b^4*c - 180
*a^4*b^2*c^2 + 125*a^5*c^3)*d^3 - 3*(9*a^3*b^5 - 51*a^4*b^3*c - 65*a^5*b*c^2)*d^2*e + 3*(3*a^4*b^4 - 22*a^5*b^
2*c - 15*a^6*c^2)*d*e^2 - (a^5*b^3 - 9*a^6*b*c)*e^3)*f)/(a^10*b^6 - 12*a^11*b^4*c + 48*a^12*b^2*c^2 - 64*a^13*
c^3)))/(a^5*b^6 - 12*a^6*b^4*c + 48*a^7*b^2*c^2 - 64*a^8*c^3))*log(-((189*b^6*c^3 - 1971*a*b^4*c^4 + 5625*a^2*
b^2*c^5 - 2500*a^3*c^6)*d^4 - (135*b^7*c^2 - 1323*a*b^5*c^3 + 2727*a^2*b^3*c^4 + 2500*a^3*b*c^5)*d^3*e + 3*(45
*a*b^6*c^2 - 558*a^2*b^4*c^3 + 1672*a^3*b^2*c^4)*d^2*e^2 - (45*a^2*b^5*c^2 - 647*a^3*b^3*c^3 + 2268*a^4*b*c^4)
*d*e^3 + (5*a^3*b^4*c^2 - 81*a^4*b^2*c^3 + 324*a^5*c^4)*e^4 - (3*a^6*b^2*c + 4*a^7*c^2)*f^4 + ((27*a^4*b^4*c +
80*a^6*c^3)*d - (9*a^5*b^3*c - 20*a^6*b*c^2)*e)*f^3 - 3*((27*a^2*b^6*c - 117*a^3*b^4*c^2 - 150*a^4*b^2*c^3 +
200*a^5*c^4)*d^2 - (18*a^3*b^5*c - 123*a^4*b^3*c^2 - 100*a^5*b*c^3)*d*e + (3*a^4*b^4*c - 28*a^5*b^2*c^2)*e^2)*
f^2 + ((81*b^8*c - 945*a*b^6*c^2 + 3213*a^2*b^4*c^3 - 3000*a^3*b^2*c^4 + 2000*a^4*c^5)*d^3 - 3*(27*a*b^7*c - 4
05*a^2*b^5*c^2 + 1461*a^3*b^3*c^3 - 500*a^4*b*c^4)*d^2*e + 3*(9*a^2*b^6*c - 165*a^3*b^4*c^2 + 692*a^4*b^2*c^3)
*d*e^2 - (3*a^3*b^5*c - 65*a^4*b^3*c^2 + 324*a^5*b*c^3)*e^3)*f)*x + 1/2*sqrt(1/2)*((27*b^11 - 486*a*b^9*c + 33
30*a^2*b^7*c^2 - 10549*a^3*b^5*c^3 + 14408*a^4*b^3*c^4 - 5200*a^5*b*c^5)*d^3 - 3*(9*a*b^10 - 177*a^2*b^8*c + 1
285*a^3*b^6*c^2 - 4138*a^4*b^4*c^3 + 5216*a^5*b^2*c^4 - 800*a^6*c^5)*d^2*e + 3*(3*a^2*b^9 - 64*a^3*b^7*c + 495
*a^4*b^5*c^2 - 1656*a^5*b^3*c^3 + 2032*a^6*b*c^4)*d*e^2 - (a^3*b^8 - 23*a^4*b^6*c + 190*a^5*b^4*c^2 - 672*a^6*
b^2*c^3 + 864*a^7*c^4)*e^3 - (a^6*b^5 - 8*a^7*b^3*c + 16*a^8*b*c^2)*f^3 + 3*((3*a^4*b^7 - 25*a^5*b^5*c + 56*a^
6*b^3*c^2 - 16*a^7*b*c^3)*d - (a^5*b^6 - 10*a^6*b^4*c + 32*a^7*b^2*c^2 - 32*a^8*c^3)*e)*f^2 - 3*((9*a^2*b^9 -
105*a^3*b^7*c + 373*a^4*b^5*c^2 - 248*a^5*b^3*c^3 - 560*a^6*b*c^4)*d^2 - 2*(3*a^3*b^8 - 40*a^4*b^6*c + 166*a^5
*b^4*c^2 - 176*a^6*b^2*c^3 - 160*a^7*c^4)*d*e + (a^4*b^7 - 15*a^5*b^5*c + 72*a^6*b^3*c^2 - 112*a^7*b*c^3)*e^2)
*f - ((3*a^5*b^10 - 55*a^6*b^8*c + 392*a^7*b^6*c^2 - 1344*a^8*b^4*c^3 + 2176*a^9*b^2*c^4 - 1280*a^10*c^5)*d -
(a^6*b^9 - 20*a^7*b^7*c + 144*a^8*b^5*c^2 - 448*a^9*b^3*c^3 + 512*a^10*b*c^4)*e - (a^7*b^8 - 8*a^8*b^6*c + 128
*a^10*b^2*c^3 - 256*a^11*c^4)*f)*sqrt((a^8*f^4 + (81*b^8 - 918*a*b^6*c + 3051*a^2*b^4*c^2 - 2550*a^3*b^2*c^3 +
625*a^4*c^4)*d^4 - 4*(27*a*b^7 - 351*a^2*b^5*c + 1197*a^3*b^3*c^2 - 550*a^4*b*c^3)*d^3*e + 6*(9*a^2*b^6 - 132
*a^3*b^4*c + 484*a^4*b^2*c^2 - 75*a^5*c^3)*d^2*e^2 - 4*(3*a^3*b^5 - 49*a^4*b^3*c + 198*a^5*b*c^2)*d*e^3 + (a^4
*b^4 - 18*a^5*b^2*c + 81*a^6*c^2)*e^4 + 4*(a^7*b*e - (3*a^6*b^2 + 5*a^7*c)*d)*f^3 + 6*((9*a^4*b^4 + 3*a^5*b^2*
c + 25*a^6*c^2)*d^2 - 2*(3*a^5*b^3 - 4*a^6*b*c)*d*e + (a^6*b^2 - 3*a^7*c)*e^2)*f^2 - 4*((27*a^2*b^6 - 108*a^3*
b^4*c - 180*a^4*b^2*c^2 + 125*a^5*c^3)*d^3 - 3*(9*a^3*b^5 - 51*a^4*b^3*c - 65*a^5*b*c^2)*d^2*e + 3*(3*a^4*b^4
- 22*a^5*b^2*c - 15*a^6*c^2)*d*e^2 - (a^5*b^3 - 9*a^6*b*c)*e^3)*f)/(a^10*b^6 - 12*a^11*b^4*c + 48*a^12*b^2*c^2
- 64*a^13*c^3)))*sqrt(-((9*b^7 - 105*a*b^5*c + 385*a^2*b^3*c^2 - 420*a^3*b*c^3)*d^2 - 2*(3*a*b^6 - 40*a^2*b^4
*c + 150*a^3*b^2*c^2 - 120*a^4*c^3)*d*e + (a^2*b^5 - 15*a^3*b^3*c + 60*a^4*b*c^2)*e^2 + (a^4*b^3 + 12*a^5*b*c)
*f^2 - 2*((3*a^2*b^5 - 13*a^3*b^3*c - 12*a^4*b*c^2)*d - (a^3*b^4 - 6*a^4*b^2*c - 24*a^5*c^2)*e)*f + (a^5*b^6 -
12*a^6*b^4*c + 48*a^7*b^2*c^2 - 64*a^8*c^3)*sqrt((a^8*f^4 + (81*b^8 - 918*a*b^6*c + 3051*a^2*b^4*c^2 - 2550*a
^3*b^2*c^3 + 625*a^4*c^4)*d^4 - 4*(27*a*b^7 - 351*a^2*b^5*c + 1197*a^3*b^3*c^2 - 550*a^4*b*c^3)*d^3*e + 6*(9*a
^2*b^6 - 132*a^3*b^4*c + 484*a^4*b^2*c^2 - 75*a^5*c^3)*d^2*e^2 - 4*(3*a^3*b^5 - 49*a^4*b^3*c + 198*a^5*b*c^2)*
d*e^3 + (a^4*b^4 - 18*a^5*b^2*c + 81*a^6*c^2)*e^4 + 4*(a^7*b*e - (3*a^6*b^2 + 5*a^7*c)*d)*f^3 + 6*((9*a^4*b^4
+ 3*a^5*b^2*c + 25*a^6*c^2)*d^2 - 2*(3*a^5*b^3 - 4*a^6*b*c)*d*e + (a^6*b^2 - 3*a^7*c)*e^2)*f^2 - 4*((27*a^2*b^
6 - 108*a^3*b^4*c - 180*a^4*b^2*c^2 + 125*a^5*c^3)*d^3 - 3*(9*a^3*b^5 - 51*a^4*b^3*c - 65*a^5*b*c^2)*d^2*e + 3
*(3*a^4*b^4 - 22*a^5*b^2*c - 15*a^6*c^2)*d*e^2 - (a^5*b^3 - 9*a^6*b*c)*e^3)*f)/(a^10*b^6 - 12*a^11*b^4*c + 48*
a^12*b^2*c^2 - 64*a^13*c^3)))/(a^5*b^6 - 12*a^6*b^4*c + 48*a^7*b^2*c^2 - 64*a^8*c^3))) - sqrt(1/2)*((a^2*b^2*c
- 4*a^3*c^2)*x^5 + (a^2*b^3 - 4*a^3*b*c)*x^3 + (a^3*b^2 - 4*a^4*c)*x)*sqrt(-((9*b^7 - 105*a*b^5*c + 385*a^2*b
^3*c^2 - 420*a^3*b*c^3)*d^2 - 2*(3*a*b^6 - 40*a^2*b^4*c + 150*a^3*b^2*c^2 - 120*a^4*c^3)*d*e + (a^2*b^5 - 15*a
^3*b^3*c + 60*a^4*b*c^2)*e^2 + (a^4*b^3 + 12*a^5*b*c)*f^2 - 2*((3*a^2*b^5 - 13*a^3*b^3*c - 12*a^4*b*c^2)*d - (
a^3*b^4 - 6*a^4*b^2*c - 24*a^5*c^2)*e)*f + (a^5*b^6 - 12*a^6*b^4*c + 48*a^7*b^2*c^2 - 64*a^8*c^3)*sqrt((a^8*f^
4 + (81*b^8 - 918*a*b^6*c + 3051*a^2*b^4*c^2 - 2550*a^3*b^2*c^3 + 625*a^4*c^4)*d^4 - 4*(27*a*b^7 - 351*a^2*b^5
*c + 1197*a^3*b^3*c^2 - 550*a^4*b*c^3)*d^3*e + 6*(9*a^2*b^6 - 132*a^3*b^4*c + 484*a^4*b^2*c^2 - 75*a^5*c^3)*d^
2*e^2 - 4*(3*a^3*b^5 - 49*a^4*b^3*c + 198*a^5*b*c^2)*d*e^3 + (a^4*b^4 - 18*a^5*b^2*c + 81*a^6*c^2)*e^4 + 4*(a^
7*b*e - (3*a^6*b^2 + 5*a^7*c)*d)*f^3 + 6*((9*a^4*b^4 + 3*a^5*b^2*c + 25*a^6*c^2)*d^2 - 2*(3*a^5*b^3 - 4*a^6*b*
c)*d*e + (a^6*b^2 - 3*a^7*c)*e^2)*f^2 - 4*((27*a^2*b^6 - 108*a^3*b^4*c - 180*a^4*b^2*c^2 + 125*a^5*c^3)*d^3 -
3*(9*a^3*b^5 - 51*a^4*b^3*c - 65*a^5*b*c^2)*d^2*e + 3*(3*a^4*b^4 - 22*a^5*b^2*c - 15*a^6*c^2)*d*e^2 - (a^5*b^3
- 9*a^6*b*c)*e^3)*f)/(a^10*b^6 - 12*a^11*b^4*c + 48*a^12*b^2*c^2 - 64*a^13*c^3)))/(a^5*b^6 - 12*a^6*b^4*c + 4
8*a^7*b^2*c^2 - 64*a^8*c^3))*log(-((189*b^6*c^3 - 1971*a*b^4*c^4 + 5625*a^2*b^2*c^5 - 2500*a^3*c^6)*d^4 - (135
*b^7*c^2 - 1323*a*b^5*c^3 + 2727*a^2*b^3*c^4 + 2500*a^3*b*c^5)*d^3*e + 3*(45*a*b^6*c^2 - 558*a^2*b^4*c^3 + 167
2*a^3*b^2*c^4)*d^2*e^2 - (45*a^2*b^5*c^2 - 647*a^3*b^3*c^3 + 2268*a^4*b*c^4)*d*e^3 + (5*a^3*b^4*c^2 - 81*a^4*b
^2*c^3 + 324*a^5*c^4)*e^4 - (3*a^6*b^2*c + 4*a^7*c^2)*f^4 + ((27*a^4*b^4*c + 80*a^6*c^3)*d - (9*a^5*b^3*c - 20
*a^6*b*c^2)*e)*f^3 - 3*((27*a^2*b^6*c - 117*a^3*b^4*c^2 - 150*a^4*b^2*c^3 + 200*a^5*c^4)*d^2 - (18*a^3*b^5*c -
123*a^4*b^3*c^2 - 100*a^5*b*c^3)*d*e + (3*a^4*b^4*c - 28*a^5*b^2*c^2)*e^2)*f^2 + ((81*b^8*c - 945*a*b^6*c^2 +
3213*a^2*b^4*c^3 - 3000*a^3*b^2*c^4 + 2000*a^4*c^5)*d^3 - 3*(27*a*b^7*c - 405*a^2*b^5*c^2 + 1461*a^3*b^3*c^3
- 500*a^4*b*c^4)*d^2*e + 3*(9*a^2*b^6*c - 165*a^3*b^4*c^2 + 692*a^4*b^2*c^3)*d*e^2 - (3*a^3*b^5*c - 65*a^4*b^3
*c^2 + 324*a^5*b*c^3)*e^3)*f)*x - 1/2*sqrt(1/2)*((27*b^11 - 486*a*b^9*c + 3330*a^2*b^7*c^2 - 10549*a^3*b^5*c^3
+ 14408*a^4*b^3*c^4 - 5200*a^5*b*c^5)*d^3 - 3*(9*a*b^10 - 177*a^2*b^8*c + 1285*a^3*b^6*c^2 - 4138*a^4*b^4*c^3
+ 5216*a^5*b^2*c^4 - 800*a^6*c^5)*d^2*e + 3*(3*a^2*b^9 - 64*a^3*b^7*c + 495*a^4*b^5*c^2 - 1656*a^5*b^3*c^3 +
2032*a^6*b*c^4)*d*e^2 - (a^3*b^8 - 23*a^4*b^6*c + 190*a^5*b^4*c^2 - 672*a^6*b^2*c^3 + 864*a^7*c^4)*e^3 - (a^6*
b^5 - 8*a^7*b^3*c + 16*a^8*b*c^2)*f^3 + 3*((3*a^4*b^7 - 25*a^5*b^5*c + 56*a^6*b^3*c^2 - 16*a^7*b*c^3)*d - (a^5
*b^6 - 10*a^6*b^4*c + 32*a^7*b^2*c^2 - 32*a^8*c^3)*e)*f^2 - 3*((9*a^2*b^9 - 105*a^3*b^7*c + 373*a^4*b^5*c^2 -
248*a^5*b^3*c^3 - 560*a^6*b*c^4)*d^2 - 2*(3*a^3*b^8 - 40*a^4*b^6*c + 166*a^5*b^4*c^2 - 176*a^6*b^2*c^3 - 160*a
^7*c^4)*d*e + (a^4*b^7 - 15*a^5*b^5*c + 72*a^6*b^3*c^2 - 112*a^7*b*c^3)*e^2)*f - ((3*a^5*b^10 - 55*a^6*b^8*c +
392*a^7*b^6*c^2 - 1344*a^8*b^4*c^3 + 2176*a^9*b^2*c^4 - 1280*a^10*c^5)*d - (a^6*b^9 - 20*a^7*b^7*c + 144*a^8*
b^5*c^2 - 448*a^9*b^3*c^3 + 512*a^10*b*c^4)*e - (a^7*b^8 - 8*a^8*b^6*c + 128*a^10*b^2*c^3 - 256*a^11*c^4)*f)*s
qrt((a^8*f^4 + (81*b^8 - 918*a*b^6*c + 3051*a^2*b^4*c^2 - 2550*a^3*b^2*c^3 + 625*a^4*c^4)*d^4 - 4*(27*a*b^7 -
351*a^2*b^5*c + 1197*a^3*b^3*c^2 - 550*a^4*b*c^3)*d^3*e + 6*(9*a^2*b^6 - 132*a^3*b^4*c + 484*a^4*b^2*c^2 - 75*
a^5*c^3)*d^2*e^2 - 4*(3*a^3*b^5 - 49*a^4*b^3*c + 198*a^5*b*c^2)*d*e^3 + (a^4*b^4 - 18*a^5*b^2*c + 81*a^6*c^2)*
e^4 + 4*(a^7*b*e - (3*a^6*b^2 + 5*a^7*c)*d)*f^3 + 6*((9*a^4*b^4 + 3*a^5*b^2*c + 25*a^6*c^2)*d^2 - 2*(3*a^5*b^3
- 4*a^6*b*c)*d*e + (a^6*b^2 - 3*a^7*c)*e^2)*f^2 - 4*((27*a^2*b^6 - 108*a^3*b^4*c - 180*a^4*b^2*c^2 + 125*a^5*
c^3)*d^3 - 3*(9*a^3*b^5 - 51*a^4*b^3*c - 65*a^5*b*c^2)*d^2*e + 3*(3*a^4*b^4 - 22*a^5*b^2*c - 15*a^6*c^2)*d*e^2
- (a^5*b^3 - 9*a^6*b*c)*e^3)*f)/(a^10*b^6 - 12*a^11*b^4*c + 48*a^12*b^2*c^2 - 64*a^13*c^3)))*sqrt(-((9*b^7 -
105*a*b^5*c + 385*a^2*b^3*c^2 - 420*a^3*b*c^3)*d^2 - 2*(3*a*b^6 - 40*a^2*b^4*c + 150*a^3*b^2*c^2 - 120*a^4*c^3
)*d*e + (a^2*b^5 - 15*a^3*b^3*c + 60*a^4*b*c^2)*e^2 + (a^4*b^3 + 12*a^5*b*c)*f^2 - 2*((3*a^2*b^5 - 13*a^3*b^3*
c - 12*a^4*b*c^2)*d - (a^3*b^4 - 6*a^4*b^2*c - 24*a^5*c^2)*e)*f + (a^5*b^6 - 12*a^6*b^4*c + 48*a^7*b^2*c^2 - 6
4*a^8*c^3)*sqrt((a^8*f^4 + (81*b^8 - 918*a*b^6*c + 3051*a^2*b^4*c^2 - 2550*a^3*b^2*c^3 + 625*a^4*c^4)*d^4 - 4*
(27*a*b^7 - 351*a^2*b^5*c + 1197*a^3*b^3*c^2 - 550*a^4*b*c^3)*d^3*e + 6*(9*a^2*b^6 - 132*a^3*b^4*c + 484*a^4*b
^2*c^2 - 75*a^5*c^3)*d^2*e^2 - 4*(3*a^3*b^5 - 49*a^4*b^3*c + 198*a^5*b*c^2)*d*e^3 + (a^4*b^4 - 18*a^5*b^2*c +
81*a^6*c^2)*e^4 + 4*(a^7*b*e - (3*a^6*b^2 + 5*a^7*c)*d)*f^3 + 6*((9*a^4*b^4 + 3*a^5*b^2*c + 25*a^6*c^2)*d^2 -
2*(3*a^5*b^3 - 4*a^6*b*c)*d*e + (a^6*b^2 - 3*a^7*c)*e^2)*f^2 - 4*((27*a^2*b^6 - 108*a^3*b^4*c - 180*a^4*b^2*c^
2 + 125*a^5*c^3)*d^3 - 3*(9*a^3*b^5 - 51*a^4*b^3*c - 65*a^5*b*c^2)*d^2*e + 3*(3*a^4*b^4 - 22*a^5*b^2*c - 15*a^
6*c^2)*d*e^2 - (a^5*b^3 - 9*a^6*b*c)*e^3)*f)/(a^10*b^6 - 12*a^11*b^4*c + 48*a^12*b^2*c^2 - 64*a^13*c^3)))/(a^5
*b^6 - 12*a^6*b^4*c + 48*a^7*b^2*c^2 - 64*a^8*c^3))) + sqrt(1/2)*((a^2*b^2*c - 4*a^3*c^2)*x^5 + (a^2*b^3 - 4*a
^3*b*c)*x^3 + (a^3*b^2 - 4*a^4*c)*x)*sqrt(-((9*b^7 - 105*a*b^5*c + 385*a^2*b^3*c^2 - 420*a^3*b*c^3)*d^2 - 2*(3
*a*b^6 - 40*a^2*b^4*c + 150*a^3*b^2*c^2 - 120*a^4*c^3)*d*e + (a^2*b^5 - 15*a^3*b^3*c + 60*a^4*b*c^2)*e^2 + (a^
4*b^3 + 12*a^5*b*c)*f^2 - 2*((3*a^2*b^5 - 13*a^3*b^3*c - 12*a^4*b*c^2)*d - (a^3*b^4 - 6*a^4*b^2*c - 24*a^5*c^2
)*e)*f - (a^5*b^6 - 12*a^6*b^4*c + 48*a^7*b^2*c^2 - 64*a^8*c^3)*sqrt((a^8*f^4 + (81*b^8 - 918*a*b^6*c + 3051*a
^2*b^4*c^2 - 2550*a^3*b^2*c^3 + 625*a^4*c^4)*d^4 - 4*(27*a*b^7 - 351*a^2*b^5*c + 1197*a^3*b^3*c^2 - 550*a^4*b*
c^3)*d^3*e + 6*(9*a^2*b^6 - 132*a^3*b^4*c + 484*a^4*b^2*c^2 - 75*a^5*c^3)*d^2*e^2 - 4*(3*a^3*b^5 - 49*a^4*b^3*
c + 198*a^5*b*c^2)*d*e^3 + (a^4*b^4 - 18*a^5*b^2*c + 81*a^6*c^2)*e^4 + 4*(a^7*b*e - (3*a^6*b^2 + 5*a^7*c)*d)*f
^3 + 6*((9*a^4*b^4 + 3*a^5*b^2*c + 25*a^6*c^2)*d^2 - 2*(3*a^5*b^3 - 4*a^6*b*c)*d*e + (a^6*b^2 - 3*a^7*c)*e^2)*
f^2 - 4*((27*a^2*b^6 - 108*a^3*b^4*c - 180*a^4*b^2*c^2 + 125*a^5*c^3)*d^3 - 3*(9*a^3*b^5 - 51*a^4*b^3*c - 65*a
^5*b*c^2)*d^2*e + 3*(3*a^4*b^4 - 22*a^5*b^2*c - 15*a^6*c^2)*d*e^2 - (a^5*b^3 - 9*a^6*b*c)*e^3)*f)/(a^10*b^6 -
12*a^11*b^4*c + 48*a^12*b^2*c^2 - 64*a^13*c^3)))/(a^5*b^6 - 12*a^6*b^4*c + 48*a^7*b^2*c^2 - 64*a^8*c^3))*log(-
((189*b^6*c^3 - 1971*a*b^4*c^4 + 5625*a^2*b^2*c^5 - 2500*a^3*c^6)*d^4 - (135*b^7*c^2 - 1323*a*b^5*c^3 + 2727*a
^2*b^3*c^4 + 2500*a^3*b*c^5)*d^3*e + 3*(45*a*b^6*c^2 - 558*a^2*b^4*c^3 + 1672*a^3*b^2*c^4)*d^2*e^2 - (45*a^2*b
^5*c^2 - 647*a^3*b^3*c^3 + 2268*a^4*b*c^4)*d*e^3 + (5*a^3*b^4*c^2 - 81*a^4*b^2*c^3 + 324*a^5*c^4)*e^4 - (3*a^6
*b^2*c + 4*a^7*c^2)*f^4 + ((27*a^4*b^4*c + 80*a^6*c^3)*d - (9*a^5*b^3*c - 20*a^6*b*c^2)*e)*f^3 - 3*((27*a^2*b^
6*c - 117*a^3*b^4*c^2 - 150*a^4*b^2*c^3 + 200*a^5*c^4)*d^2 - (18*a^3*b^5*c - 123*a^4*b^3*c^2 - 100*a^5*b*c^3)*
d*e + (3*a^4*b^4*c - 28*a^5*b^2*c^2)*e^2)*f^2 + ((81*b^8*c - 945*a*b^6*c^2 + 3213*a^2*b^4*c^3 - 3000*a^3*b^2*c
^4 + 2000*a^4*c^5)*d^3 - 3*(27*a*b^7*c - 405*a^2*b^5*c^2 + 1461*a^3*b^3*c^3 - 500*a^4*b*c^4)*d^2*e + 3*(9*a^2*
b^6*c - 165*a^3*b^4*c^2 + 692*a^4*b^2*c^3)*d*e^2 - (3*a^3*b^5*c - 65*a^4*b^3*c^2 + 324*a^5*b*c^3)*e^3)*f)*x +
1/2*sqrt(1/2)*((27*b^11 - 486*a*b^9*c + 3330*a^2*b^7*c^2 - 10549*a^3*b^5*c^3 + 14408*a^4*b^3*c^4 - 5200*a^5*b*
c^5)*d^3 - 3*(9*a*b^10 - 177*a^2*b^8*c + 1285*a^3*b^6*c^2 - 4138*a^4*b^4*c^3 + 5216*a^5*b^2*c^4 - 800*a^6*c^5)
*d^2*e + 3*(3*a^2*b^9 - 64*a^3*b^7*c + 495*a^4*b^5*c^2 - 1656*a^5*b^3*c^3 + 2032*a^6*b*c^4)*d*e^2 - (a^3*b^8 -
23*a^4*b^6*c + 190*a^5*b^4*c^2 - 672*a^6*b^2*c^3 + 864*a^7*c^4)*e^3 - (a^6*b^5 - 8*a^7*b^3*c + 16*a^8*b*c^2)*
f^3 + 3*((3*a^4*b^7 - 25*a^5*b^5*c + 56*a^6*b^3*c^2 - 16*a^7*b*c^3)*d - (a^5*b^6 - 10*a^6*b^4*c + 32*a^7*b^2*c
^2 - 32*a^8*c^3)*e)*f^2 - 3*((9*a^2*b^9 - 105*a^3*b^7*c + 373*a^4*b^5*c^2 - 248*a^5*b^3*c^3 - 560*a^6*b*c^4)*d
^2 - 2*(3*a^3*b^8 - 40*a^4*b^6*c + 166*a^5*b^4*c^2 - 176*a^6*b^2*c^3 - 160*a^7*c^4)*d*e + (a^4*b^7 - 15*a^5*b^
5*c + 72*a^6*b^3*c^2 - 112*a^7*b*c^3)*e^2)*f + ((3*a^5*b^10 - 55*a^6*b^8*c + 392*a^7*b^6*c^2 - 1344*a^8*b^4*c^
3 + 2176*a^9*b^2*c^4 - 1280*a^10*c^5)*d - (a^6*b^9 - 20*a^7*b^7*c + 144*a^8*b^5*c^2 - 448*a^9*b^3*c^3 + 512*a^
10*b*c^4)*e - (a^7*b^8 - 8*a^8*b^6*c + 128*a^10*b^2*c^3 - 256*a^11*c^4)*f)*sqrt((a^8*f^4 + (81*b^8 - 918*a*b^6
*c + 3051*a^2*b^4*c^2 - 2550*a^3*b^2*c^3 + 625*a^4*c^4)*d^4 - 4*(27*a*b^7 - 351*a^2*b^5*c + 1197*a^3*b^3*c^2 -
550*a^4*b*c^3)*d^3*e + 6*(9*a^2*b^6 - 132*a^3*b^4*c + 484*a^4*b^2*c^2 - 75*a^5*c^3)*d^2*e^2 - 4*(3*a^3*b^5 -
49*a^4*b^3*c + 198*a^5*b*c^2)*d*e^3 + (a^4*b^4 - 18*a^5*b^2*c + 81*a^6*c^2)*e^4 + 4*(a^7*b*e - (3*a^6*b^2 + 5*
a^7*c)*d)*f^3 + 6*((9*a^4*b^4 + 3*a^5*b^2*c + 25*a^6*c^2)*d^2 - 2*(3*a^5*b^3 - 4*a^6*b*c)*d*e + (a^6*b^2 - 3*a
^7*c)*e^2)*f^2 - 4*((27*a^2*b^6 - 108*a^3*b^4*c - 180*a^4*b^2*c^2 + 125*a^5*c^3)*d^3 - 3*(9*a^3*b^5 - 51*a^4*b
^3*c - 65*a^5*b*c^2)*d^2*e + 3*(3*a^4*b^4 - 22*a^5*b^2*c - 15*a^6*c^2)*d*e^2 - (a^5*b^3 - 9*a^6*b*c)*e^3)*f)/(
a^10*b^6 - 12*a^11*b^4*c + 48*a^12*b^2*c^2 - 64*a^13*c^3)))*sqrt(-((9*b^7 - 105*a*b^5*c + 385*a^2*b^3*c^2 - 42
0*a^3*b*c^3)*d^2 - 2*(3*a*b^6 - 40*a^2*b^4*c + 150*a^3*b^2*c^2 - 120*a^4*c^3)*d*e + (a^2*b^5 - 15*a^3*b^3*c +
60*a^4*b*c^2)*e^2 + (a^4*b^3 + 12*a^5*b*c)*f^2 - 2*((3*a^2*b^5 - 13*a^3*b^3*c - 12*a^4*b*c^2)*d - (a^3*b^4 - 6
*a^4*b^2*c - 24*a^5*c^2)*e)*f - (a^5*b^6 - 12*a^6*b^4*c + 48*a^7*b^2*c^2 - 64*a^8*c^3)*sqrt((a^8*f^4 + (81*b^8
- 918*a*b^6*c + 3051*a^2*b^4*c^2 - 2550*a^3*b^2*c^3 + 625*a^4*c^4)*d^4 - 4*(27*a*b^7 - 351*a^2*b^5*c + 1197*a
^3*b^3*c^2 - 550*a^4*b*c^3)*d^3*e + 6*(9*a^2*b^6 - 132*a^3*b^4*c + 484*a^4*b^2*c^2 - 75*a^5*c^3)*d^2*e^2 - 4*(
3*a^3*b^5 - 49*a^4*b^3*c + 198*a^5*b*c^2)*d*e^3 + (a^4*b^4 - 18*a^5*b^2*c + 81*a^6*c^2)*e^4 + 4*(a^7*b*e - (3*
a^6*b^2 + 5*a^7*c)*d)*f^3 + 6*((9*a^4*b^4 + 3*a^5*b^2*c + 25*a^6*c^2)*d^2 - 2*(3*a^5*b^3 - 4*a^6*b*c)*d*e + (a
^6*b^2 - 3*a^7*c)*e^2)*f^2 - 4*((27*a^2*b^6 - 108*a^3*b^4*c - 180*a^4*b^2*c^2 + 125*a^5*c^3)*d^3 - 3*(9*a^3*b^
5 - 51*a^4*b^3*c - 65*a^5*b*c^2)*d^2*e + 3*(3*a^4*b^4 - 22*a^5*b^2*c - 15*a^6*c^2)*d*e^2 - (a^5*b^3 - 9*a^6*b*
c)*e^3)*f)/(a^10*b^6 - 12*a^11*b^4*c + 48*a^12*b^2*c^2 - 64*a^13*c^3)))/(a^5*b^6 - 12*a^6*b^4*c + 48*a^7*b^2*c
^2 - 64*a^8*c^3))) - sqrt(1/2)*((a^2*b^2*c - 4*a^3*c^2)*x^5 + (a^2*b^3 - 4*a^3*b*c)*x^3 + (a^3*b^2 - 4*a^4*c)*
x)*sqrt(-((9*b^7 - 105*a*b^5*c + 385*a^2*b^3*c^2 - 420*a^3*b*c^3)*d^2 - 2*(3*a*b^6 - 40*a^2*b^4*c + 150*a^3*b^
2*c^2 - 120*a^4*c^3)*d*e + (a^2*b^5 - 15*a^3*b^3*c + 60*a^4*b*c^2)*e^2 + (a^4*b^3 + 12*a^5*b*c)*f^2 - 2*((3*a^
2*b^5 - 13*a^3*b^3*c - 12*a^4*b*c^2)*d - (a^3*b^4 - 6*a^4*b^2*c - 24*a^5*c^2)*e)*f - (a^5*b^6 - 12*a^6*b^4*c +
48*a^7*b^2*c^2 - 64*a^8*c^3)*sqrt((a^8*f^4 + (81*b^8 - 918*a*b^6*c + 3051*a^2*b^4*c^2 - 2550*a^3*b^2*c^3 + 62
5*a^4*c^4)*d^4 - 4*(27*a*b^7 - 351*a^2*b^5*c + 1197*a^3*b^3*c^2 - 550*a^4*b*c^3)*d^3*e + 6*(9*a^2*b^6 - 132*a^
3*b^4*c + 484*a^4*b^2*c^2 - 75*a^5*c^3)*d^2*e^2 - 4*(3*a^3*b^5 - 49*a^4*b^3*c + 198*a^5*b*c^2)*d*e^3 + (a^4*b^
4 - 18*a^5*b^2*c + 81*a^6*c^2)*e^4 + 4*(a^7*b*e - (3*a^6*b^2 + 5*a^7*c)*d)*f^3 + 6*((9*a^4*b^4 + 3*a^5*b^2*c +
25*a^6*c^2)*d^2 - 2*(3*a^5*b^3 - 4*a^6*b*c)*d*e + (a^6*b^2 - 3*a^7*c)*e^2)*f^2 - 4*((27*a^2*b^6 - 108*a^3*b^4
*c - 180*a^4*b^2*c^2 + 125*a^5*c^3)*d^3 - 3*(9*a^3*b^5 - 51*a^4*b^3*c - 65*a^5*b*c^2)*d^2*e + 3*(3*a^4*b^4 - 2
2*a^5*b^2*c - 15*a^6*c^2)*d*e^2 - (a^5*b^3 - 9*a^6*b*c)*e^3)*f)/(a^10*b^6 - 12*a^11*b^4*c + 48*a^12*b^2*c^2 -
64*a^13*c^3)))/(a^5*b^6 - 12*a^6*b^4*c + 48*a^7*b^2*c^2 - 64*a^8*c^3))*log(-((189*b^6*c^3 - 1971*a*b^4*c^4 + 5
625*a^2*b^2*c^5 - 2500*a^3*c^6)*d^4 - (135*b^7*c^2 - 1323*a*b^5*c^3 + 2727*a^2*b^3*c^4 + 2500*a^3*b*c^5)*d^3*e
+ 3*(45*a*b^6*c^2 - 558*a^2*b^4*c^3 + 1672*a^3*b^2*c^4)*d^2*e^2 - (45*a^2*b^5*c^2 - 647*a^3*b^3*c^3 + 2268*a^
4*b*c^4)*d*e^3 + (5*a^3*b^4*c^2 - 81*a^4*b^2*c^3 + 324*a^5*c^4)*e^4 - (3*a^6*b^2*c + 4*a^7*c^2)*f^4 + ((27*a^4
*b^4*c + 80*a^6*c^3)*d - (9*a^5*b^3*c - 20*a^6*b*c^2)*e)*f^3 - 3*((27*a^2*b^6*c - 117*a^3*b^4*c^2 - 150*a^4*b^
2*c^3 + 200*a^5*c^4)*d^2 - (18*a^3*b^5*c - 123*a^4*b^3*c^2 - 100*a^5*b*c^3)*d*e + (3*a^4*b^4*c - 28*a^5*b^2*c^
2)*e^2)*f^2 + ((81*b^8*c - 945*a*b^6*c^2 + 3213*a^2*b^4*c^3 - 3000*a^3*b^2*c^4 + 2000*a^4*c^5)*d^3 - 3*(27*a*b
^7*c - 405*a^2*b^5*c^2 + 1461*a^3*b^3*c^3 - 500*a^4*b*c^4)*d^2*e + 3*(9*a^2*b^6*c - 165*a^3*b^4*c^2 + 692*a^4*
b^2*c^3)*d*e^2 - (3*a^3*b^5*c - 65*a^4*b^3*c^2 + 324*a^5*b*c^3)*e^3)*f)*x - 1/2*sqrt(1/2)*((27*b^11 - 486*a*b^
9*c + 3330*a^2*b^7*c^2 - 10549*a^3*b^5*c^3 + 14408*a^4*b^3*c^4 - 5200*a^5*b*c^5)*d^3 - 3*(9*a*b^10 - 177*a^2*b
^8*c + 1285*a^3*b^6*c^2 - 4138*a^4*b^4*c^3 + 5216*a^5*b^2*c^4 - 800*a^6*c^5)*d^2*e + 3*(3*a^2*b^9 - 64*a^3*b^7
*c + 495*a^4*b^5*c^2 - 1656*a^5*b^3*c^3 + 2032*a^6*b*c^4)*d*e^2 - (a^3*b^8 - 23*a^4*b^6*c + 190*a^5*b^4*c^2 -
672*a^6*b^2*c^3 + 864*a^7*c^4)*e^3 - (a^6*b^5 - 8*a^7*b^3*c + 16*a^8*b*c^2)*f^3 + 3*((3*a^4*b^7 - 25*a^5*b^5*c
+ 56*a^6*b^3*c^2 - 16*a^7*b*c^3)*d - (a^5*b^6 - 10*a^6*b^4*c + 32*a^7*b^2*c^2 - 32*a^8*c^3)*e)*f^2 - 3*((9*a^
2*b^9 - 105*a^3*b^7*c + 373*a^4*b^5*c^2 - 248*a^5*b^3*c^3 - 560*a^6*b*c^4)*d^2 - 2*(3*a^3*b^8 - 40*a^4*b^6*c +
166*a^5*b^4*c^2 - 176*a^6*b^2*c^3 - 160*a^7*c^4)*d*e + (a^4*b^7 - 15*a^5*b^5*c + 72*a^6*b^3*c^2 - 112*a^7*b*c
^3)*e^2)*f + ((3*a^5*b^10 - 55*a^6*b^8*c + 392*a^7*b^6*c^2 - 1344*a^8*b^4*c^3 + 2176*a^9*b^2*c^4 - 1280*a^10*c
^5)*d - (a^6*b^9 - 20*a^7*b^7*c + 144*a^8*b^5*c^2 - 448*a^9*b^3*c^3 + 512*a^10*b*c^4)*e - (a^7*b^8 - 8*a^8*b^6
*c + 128*a^10*b^2*c^3 - 256*a^11*c^4)*f)*sqrt((a^8*f^4 + (81*b^8 - 918*a*b^6*c + 3051*a^2*b^4*c^2 - 2550*a^3*b
^2*c^3 + 625*a^4*c^4)*d^4 - 4*(27*a*b^7 - 351*a^2*b^5*c + 1197*a^3*b^3*c^2 - 550*a^4*b*c^3)*d^3*e + 6*(9*a^2*b
^6 - 132*a^3*b^4*c + 484*a^4*b^2*c^2 - 75*a^5*c^3)*d^2*e^2 - 4*(3*a^3*b^5 - 49*a^4*b^3*c + 198*a^5*b*c^2)*d*e^
3 + (a^4*b^4 - 18*a^5*b^2*c + 81*a^6*c^2)*e^4 + 4*(a^7*b*e - (3*a^6*b^2 + 5*a^7*c)*d)*f^3 + 6*((9*a^4*b^4 + 3*
a^5*b^2*c + 25*a^6*c^2)*d^2 - 2*(3*a^5*b^3 - 4*a^6*b*c)*d*e + (a^6*b^2 - 3*a^7*c)*e^2)*f^2 - 4*((27*a^2*b^6 -
108*a^3*b^4*c - 180*a^4*b^2*c^2 + 125*a^5*c^3)*d^3 - 3*(9*a^3*b^5 - 51*a^4*b^3*c - 65*a^5*b*c^2)*d^2*e + 3*(3*
a^4*b^4 - 22*a^5*b^2*c - 15*a^6*c^2)*d*e^2 - (a^5*b^3 - 9*a^6*b*c)*e^3)*f)/(a^10*b^6 - 12*a^11*b^4*c + 48*a^12
*b^2*c^2 - 64*a^13*c^3)))*sqrt(-((9*b^7 - 105*a*b^5*c + 385*a^2*b^3*c^2 - 420*a^3*b*c^3)*d^2 - 2*(3*a*b^6 - 40
*a^2*b^4*c + 150*a^3*b^2*c^2 - 120*a^4*c^3)*d*e + (a^2*b^5 - 15*a^3*b^3*c + 60*a^4*b*c^2)*e^2 + (a^4*b^3 + 12*
a^5*b*c)*f^2 - 2*((3*a^2*b^5 - 13*a^3*b^3*c - 12*a^4*b*c^2)*d - (a^3*b^4 - 6*a^4*b^2*c - 24*a^5*c^2)*e)*f - (a
^5*b^6 - 12*a^6*b^4*c + 48*a^7*b^2*c^2 - 64*a^8*c^3)*sqrt((a^8*f^4 + (81*b^8 - 918*a*b^6*c + 3051*a^2*b^4*c^2
- 2550*a^3*b^2*c^3 + 625*a^4*c^4)*d^4 - 4*(27*a*b^7 - 351*a^2*b^5*c + 1197*a^3*b^3*c^2 - 550*a^4*b*c^3)*d^3*e
+ 6*(9*a^2*b^6 - 132*a^3*b^4*c + 484*a^4*b^2*c^2 - 75*a^5*c^3)*d^2*e^2 - 4*(3*a^3*b^5 - 49*a^4*b^3*c + 198*a^5
*b*c^2)*d*e^3 + (a^4*b^4 - 18*a^5*b^2*c + 81*a^6*c^2)*e^4 + 4*(a^7*b*e - (3*a^6*b^2 + 5*a^7*c)*d)*f^3 + 6*((9*
a^4*b^4 + 3*a^5*b^2*c + 25*a^6*c^2)*d^2 - 2*(3*a^5*b^3 - 4*a^6*b*c)*d*e + (a^6*b^2 - 3*a^7*c)*e^2)*f^2 - 4*((2
7*a^2*b^6 - 108*a^3*b^4*c - 180*a^4*b^2*c^2 + 125*a^5*c^3)*d^3 - 3*(9*a^3*b^5 - 51*a^4*b^3*c - 65*a^5*b*c^2)*d
^2*e + 3*(3*a^4*b^4 - 22*a^5*b^2*c - 15*a^6*c^2)*d*e^2 - (a^5*b^3 - 9*a^6*b*c)*e^3)*f)/(a^10*b^6 - 12*a^11*b^4
*c + 48*a^12*b^2*c^2 - 64*a^13*c^3)))/(a^5*b^6 - 12*a^6*b^4*c + 48*a^7*b^2*c^2 - 64*a^8*c^3))) - 4*(a*b^2 - 4*
a^2*c)*d)/((a^2*b^2*c - 4*a^3*c^2)*x^5 + (a^2*b^3 - 4*a^3*b*c)*x^3 + (a^3*b^2 - 4*a^4*c)*x)

________________________________________________________________________________________

Sympy [F(-1)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Timed out} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((f*x**4+e*x**2+d)/x**2/(c*x**4+b*x**2+a)**2,x)

[Out]

Timed out

________________________________________________________________________________________

Giac [F(-2)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Exception raised: NotImplementedError} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((f*x^4+e*x^2+d)/x^2/(c*x^4+b*x^2+a)^2,x, algorithm="giac")

[Out]

Exception raised: NotImplementedError